Criticality in one dimension with inverse square-law potentials 
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It is demonstrated that the scaled order parameter for ferromagnetic Ising and three-state Potts 
chains with inverse-square interactions exhibits a universal critical jump, in analogy with the su- 
perfluid density in helium films. Renormalization-group arguments are combined with numerical 
simulations ol systems containing up to one million lattice sites to accurately determine the critical 
properties ol these models. In strong contrast with earlier work, compelling quantitative evidence 
for the Kosterlitz-Thouless-like character ol the phase transition is provided. 
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One-dimensional scalar models with inverse-square po- 
tentials have held a special position in the field of statisti- 
cal mechanics for several decades. The surprisingly rich 
behavior emerging from these conceptually simple sys- 
tems has attracted the attention of a variety of workers 
in the field. The work of Ruelle ]l| and Dyson Q proved 
that pair interactions decaying as 1 jr 2 (where r denotes 
the distance between two local variables) represent the 
boundary between models with and without a phase tran- 
sition at nonzero temperature; but it took until the 1980s 
before Frohlich and Spencer |}| proved rigorously the ex- 
istence of a phase transition for the borderline case it- 
self. Meanwhile, early work by Anderson and Yuval Q 
(who were motivated by the connection to the Kondo 
problem) had suggested that this phase transition has a 
special character, with an order parameter that exhibits 
both a discontinuity and a superposed singularity at the 
critical point. This behavior later gained widespread at- 
tention through the Kosterlitz-Thouless (KT) transition 
in the two-dimensional (2D) XY model fa]. Indeed, the 
vortex-unbinding scenario in the latter case has a direct 
analog in terms of logarithmically interacting defects (do- 
main walls) in the "kink-gas" representation of the Ising 
chain considered in [EJ . This treatment in the spirit of the 
renormalization-group (RG) theory was extended in || 
and most notably by Cardy Q], who demonstrated that 
it can be generalized to large classes of discrete models. 
In this sense, these one-dimensional models can be con- 
sidered as the simplest examples of the class of "topo- 
logical" phase transitions, which comprises phenomena 
like 2D superfluidity and superconductivity, as well as 
dislocation-mediated melting of crystals. 

An exact lower bound on the critical order-parameter 
jump has been obtained by Aizenman et ol. || for Ising 
and Potts chains with 1/r 2 interactions, and the existence 



of an intermediate ordered phase, where the correlation 
function decays like a temperature-dependent power-law, 
has been proven by Imbrie and Newman Despite 
these rigorous results the literature on these models con- 
tains various contradictory claims and unconfirmed con- 
jectures, while few of the RG-based predictions for the 
critical properties have been verified by alternative meth- 
ods: Thus (i) the critical temperature T c has been esti- 
mated by a great variety of techniques, but little consis- 
tency has been obtained until now, with estimates dif- 
fering by up to 30%. Indeed, there is even a conjec- 
ture for T c of the g-state Potts chain that is independent 
of q ]io| ], despite the fact that this appears to contradict 
intuitive energy-entropy arguments, (ii) In contrast, it 
has been claimed on the basis of numerical simulations 
that the transition in the g-state Potts chain is contin- 
uous for q < q c and of first order for q > q c [fj"l| , al- 
though this in turn disagrees with the KT-like transition 
(with g-dependent exponents) found by Cardy @ . (iii) It 
has been conjectured that the power-law exponent of the 
correlation function in the intermediate ordered phase is 
directly related to the value of the order parameter [||. 
(iv) Numerical simulations have been unable to conclude 
whether the specific heat diverges at T c or even whether 
it displays a maximum at T c Jl2| . (v) No numerical work 
for a one-dimensional system has been able to reproduce 
one of the actual hallmarks of the KT universality class, 
namely the predicted exponential divergence of the cor- 
relation length and the susceptibility for T \T C . 

It is the aim of the present work to settle these open 
questions in a convincing way, although (ii) will be ad- 
dressed only partially, by combining finite-size scaling re- 
lations obtained from the RG predictions (which until 
now have been ignored in nearly all numerical studies) 
with accurate numerical simulations that span an un- 
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precedented range in system sizes. Indeed, the long-range 
nature of the interactions has posed a major hurdle for 
the numerical study of these models: Whereas almost 
two decades ago Monte Carlo (MC) simulations could 
reach a system of L = 256 lattice sites [Q, the most re- 
cent work still could only attain L = 900 O]. By using 
appropriate techniques, we have now been able to study 
systems as large as L = 10 6 . This allows us not only 
to answer questions intrinsic to the one-dimensional sys- 
tems under study, but also to obtain what we believe are 
the most accurate numerical results to date for a system 
exhibiting a KT-like phase transition: Indeed, although 
the 2D XY and Villain models have been examined in 
many numerical studies, it has proven extremely difficult 
to obtain full consistency with the RG predictions [^3| . 
Because of the exponential divergence of the correlation 
length £, it is essential to cover a large range in linear sys- 
tem sizes, which makes a one-dimensional system easier 
to study than a 2D one. In our analysis, we start from the 
hypothesis that the structure of the scaling functions is 
correctly described by RG theory, and examine the con- 
sistency between our numerical results and the scaling 
predictions to test the correctness of this hypothesis. 

We have first focused attention on the Ising chain (or 
equivalently the q = 2 Potts chain), described by the re- 
duced Hamiltonian f¥H = —KY^uj\ s i s j/ r iji where the 
sum runs over all spin pairs, s,; = ±1, and = — r^j . 
We have carried out accurate MC simulations for chains 
of L lattice sites, with 10 < L < 10 6 , over a large tem- 
perature range. These results could be obtained owing 
to the use of a cluster algorithm that suppresses critical 
slowing down and has an efficiency that is independent 
of the interaction range |14|. 



Although the critical coupling K c is a nonuniversal 
quantity, the large number of available estimates makes 
it interesting in its own right. In addition, its precise 
value proves indispensable if one wants to address the 
issues mentioned above. Since finite-size corrections de- 
cay only logarithmically near a KT transition, the de- 
termination of K c is notoriously difficult, and a new ap- 
proach had been found. To this end, we have relied on 
the observation that the properly scaled order param- 
eter 'J = K(m 2 ), with m = L^ 1 J2f =1 s^, is expected 
to exhibit a universal jump at the critical point. In- 
deed, as mentioned, it has been shown rigorously J^] that 
^{K > K c ) > \ [whereas V(K < K c ) = 0]. Further- 
more, it was conjectured in Q that this inequality is sat- 
urated at K = K c (which seems to be implied by the 
arguments of [Q, but note the comments in ||). While 
it can be easily shown that 4" is invariant under a (real- 
space) RG transformation of the kink-gas representation 
of the Ising chain [ jl5| , it appears not to have been ex- 
plicitly realized that the critical jump ^ c = ^{K c ) must 
be universal, in the sense that it will have the same mag- 
nitude in every system with couplings that decay asymp- 



totically like r~ 2 . We remark that also for the analogous 
phenomenon in two dimensions (the celebrated universal 
jump in the superfluid density of a 4 He film) this uni- 



versality was only realized 
by Kosterlitz and Thouless [5] 
discontinuity, <3>(if > K, 
perposed singularity: 



after the original work 
In addition to the jump 
is predicted to exhibit a su- 



V(K) = ^ c + ai (K- K c f + a 2 {K- K c 



, (1) 



with v = i for the Ising chain; the ellipsis also includes 
nonsingular correction terms. In what follows, we demon- 
strate how to exploit the scaling properties of ^> to accu- 
rately locate K c . Figure g shows ^(K,L) as a function 
of K for different system sizes. While 'J is already sharp- 
ening up for relatively small L, it is only for very large 
systems that one observes the development of a discon- 
tinuity and an additional singularity. 

To proceed further, we recall the lowest-order RG flow 
equations j|,D for the coupling and the defect fugacity y, 
namely 



dK/dl = -4K y 2 , 
dy/dl = -y(2K - 1) , 



(2) 
(3) 



which (under the additional assumption that the system 
size scales with the rescaling parameter e l ) allow us to 
derive the finite-size scaling behavior of ^>(K,L) in the 
low-temperature regime. Since this derivation is some- 
what involved |l5j| , we only quote the final result here, 
namely 



V(K,L) = *(if,oo) 



1 + hL- 2 ^ 



(4) 



where = ^(K, oo)/ 1 $ c . Nonsingular corrections be- 
come important only when ^ > 2^> c . The appear- 
ance of ^{K, oo) both as an amplitude and as an ex- 
ponent is typical for a KT transition. A least-squares 
fit of this equation to ^(K,L) has been used to ob- 
tain ^(K, oo) for several values of K, as shown in the 
inset in Fig. @. This also yielded * c = 0.49 ± 0.02. 
Keeping ^f c fixed at the predicted value | |§[lf| yielded 
K c = 0.6552 (2). Subsequently, Eq. (0) allowed us to 
determine v = 0.52 ± 0.03. We have supplemented this 
low-temperature approach with an analysis of data ly- 
ing within the high-temperature finite-size regime, where 
L <C C an d K < K c . The algebraic corrections in (TJ) now 
become powers of 1/ In L and temperature-dependent 
terms take the form L/£, with £ = exp[S / \K-K c f] jij§. 
We then find \& c = 0.496 (3), in excellent agreement with 
the RG prediction, and K c = 0.6548 (14). Keeping * c 
now fixed at \ yields K c = 0.6555 (4), very close to the 
value resulting from the low-temperature approach. 

Independent estimates for K c and v can be obtained 
from the high-temperature susceptibility, which should 
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display the same exponential divergence as the correla- 
tion length. RG considerations regarding higher-order 
contributions |jl7) suggest x ~ ^p[B Q (K — K c )~ u + 
B 2 (K - K c y + 0{{K - K c ) 29 )}. Numerical results for 
L = 4 • 10 5 and L = 10 6 are shown in Fig. [| along 
with a truncated high-tcmperature expansion of order 
0(K 8 ) and a fit to the RG expression. The agree- 
ment with the series expansion is good, but breaks down 
for K > 0.3. The RG expression provides an excellent 
description of the MC data for 0.2 < K < 0.61, con- 
firming the predicted exponential divergence over a large 
temperature range. Owing to the covariance of the fit pa- 
rameters, it is difficult to determine them simultaneously: 
We obtained v = 0.54 ± 0.03 as our best estimate. How- 
ever, keeping v = \ leads to K c — 0.6551 (6), which we 
view as a strong confirmation of our analysis of $ . The 
rounding of x close to K c (see Fig. ||) is due to the finite 
system size. We estimate that the correlation length £ is 
> 10 here, which may be compared with a maximum of 
£ w 140 reached in recent d — 2 simulations jl3| . 

Subsequently, we carried out a similar analysis of 
VP and x f° r the g-state Potts model with q = 3, 
which has a Hamiltonian of the form /37i Potts = 
-K Potts <W4> with Si = l,...,q. Cardy § has 
found that, for positive integer q, this system can be de- 
scribed by RG equations that are very similar to those 
for the q = 2 (Ising) case, but yielding v — 2/(q + 2). 
This implies that the entire analysis of the critical be- 
havior can be carried through as before |l5): ^ should 
still exhibit a universal jump (of magnitude 1 in units 
of K Potts ) with a superposed singularity as in (|l|) and x 
should again display an exponential singularity. The or- 
der parameter used for the determination of W and x. is 
given by m Potts = max^!,...,^ I]f =1 (5 SjCe - l)]/(<7 - 1), 
with q > 1. An analysis of \1/ for system sizes up to 
L = 10 6 yielded K c = 1.4105 (2), * c = 1.02 (2), and 
v(q = 3) = 0.41 (2). A least-squares fit of the data 
for x was consistent with an exponential divergence with 
v(q = 3) = 0.42 (2). Not only is ty c in good agreement 
with the predicted universal value, but to our knowl- 
edge this is also the very first independent corrobora- 
tion of the generalized RG scenario of Ref. §, which 
predicts that v(q = 3) = |. As expected from en- 
tropic considerations, K Potts (q = 3) is somewhat higher 
than K Potts (q = 2) = 1.3104 (4), which in turn re- 
futes two conjectures based on a different real-space RG 
analysis jnj, namely that (i) K Potts would be indepen- 
dent of q and (ii) A"^ otts would be exactly given by 
12/7T 2 = 1.21585 . . . We note that our results are in strong 
contrast with earlier numerical work: Both a transfer- 
matrix study and a very recent MC study for q — 2 |h| 
and q = 3 [H||^] concluded that the correlation length 
diverges algebraically (rather than exponentially) with a 
q-dependent power-law exponent. 

We now return to some interesting aspects of the Ising 
chain. First, it has been concluded in Ref. Q that 



the critical spin-spin correlation function should decay 
logarithmically, g(r) = (sQS r ) — (sq) 2 ~ 1/lnr, where 
now (sq) 7^ 0, unlike the situation for a typical criti- 
cal phase transition. As illustrated in Fig. ||, we have 
been able to corroborate this behavior over a consid- 
erable range in r. In the high-temperature phase, the 
decay of g(r) is bounded by the spin-spin interaction, 
but for K > K c it is expected that g(r) ~ r~ e , with 
9 = AyJ{K - K c )/K c Q. We note that this behavior 
would be consistent with the conjecture of Ref. viz. 
that 8 = min[4(* - ±),2] for K > K c . We have de- 
termined 9 for a large number of couplings, see Fig. |^, 
and have subsequently fitted these to a power series in 
^JK - K c , finding 9 = (1.81±0M)^(K - K c )/K c . This 
agrees with Aai\fW c = 1.73 (4), where a\ is defined via 
Eq. (Q), corroborating the conjectured relation with "P. 
However, the prefactor is approximately half as large as 
predicted in Ref. |l2| . Farther away from K c it becomes 
increasingly difficult to obtain accurate results for 9(K), 
owing to the rapid decay of g{r). However, since our re- 
sults corroborate the presumed relation between 9 and VP, 
we exploit this relation to determine two special values 
of K, namely the coupling strength K\ for which 9 = 1 
and the smallest coupling K2 > K c for which 9 = 2. K\ 
has the special property that the magnetic susceptibility 
will diverge over the full coupling range K c < K < K±, 
as follows from the fluctuation-dissipation theorem (and 
as is actually confirmed by the numerical results |p5[). 
We find Kx = 0.78 (1) and K 2 = 1.02 (3). Thus, the RG 
prediction K\ = j^K c « 0.696 [[l2| is too low, due to the 
difference in the prefactor for 9(K). K 2 marks the end 
of the regime where g(r) displays a leading temperature- 
dependent exponent. 

Finally, we consider the specific heat C, which is cal- 
culated from the fluctuations in the internal energy. The 
efficiency of the original MC algorithm JlJ] is greatly re- 
duced when this energy has to be calculated (although 
recently this problem has largely been overcome | f20[ ) 
and hence only much smaller systems could be accessed. 
Nevertheless, this turns out not to be a limiting factor, 
since C has already converged within the range of ac- 
cessible system sizes: Fig. ^ shows unambiguously that 
C is a nondivergent quantity which takes its maximum 
at a temperature that lies approximately 14% above the 
critical temperature. 

In summary, we have provided compelling evidence 
that the inverse-square Ising and q = 3 Potts chain 
display the one-dimensional analog of a Kosterlitz- 
Thouless transition, in disagreement with earlier numer- 
ical studies, but in accurate quantitative agreement with 
renormalization-group predictions. The (exponentially 
diverging) correlation length could be examined over 
more than five orders of magnitude, making this the most 
tractable KT-like transition. 

It is a great pleasure to acknowledge stimulating re- 
marks by Prof. Kurt Binder and helpful comments by 



3 



Prof. Michael E. Fisher. E.L. gratefully acknowledges 
the hospitality of the Condensed Matter Theory Group at 
the Johannes Gutenberg-Universitat Mainz, where most 
of this work has been carried out, and support from the 
National Science Foundation (through Grant No. CHE 
99-81772 to M.E. Fisher) and from the Department of 
Energy, Office of Basic Energy Sciences (through Grant 
No. DE-FG02-98ER14858 to A.Z. Panagiotopoulos). 



[1] D. Ruelle, Comm. Math. Phys. 9, 267 (1968). 
[2] F. J. Dyson, Comm. Math. Phys. 12, 91 (1969). 
[3] J. Frohlich and T. Spencer, Comm. Math. Phys. 84, 87 
(1982). 

[4] P. W. Anderson and G. Yuval, J. Phys. C 4, 607 (1971). 

[5] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 
(1973); J. M. Kosterlitz, J. Phys. C 7, 1046 (1974). 

[6] J. M. Kosterlitz, Phys. Rev. Lett. 37, 1577 (1976). 

[7] J. L. Cardy, J. Phys. A 14, 1407 (1981). 

[8] M. Aizenman, J. T. Chayes, L. Chayes, and C. M. New- 
man, J. Stat. Phys. 50, 1 (1988). Note that this paper is 
concerned with the Potts chain, and hence the coupling 
strengths are multiplied by a factor 2, K Potts = 2K Ia[as . 

[9] J. Z. Imbrie and C. M. Newman, Comm. Math. Phys. 
118, 303 (1988). 
[10] S. A. Cannas and A. C. N. de Magalhaes, J. Phys. A 30, 
3345 (1997). 

[11] E. Bayong, H. T. Diep, and V. Dotsenko, Phys. Rev. 

Lett. 83, 14 (1999). 
[12] J. Bhattacharjee, S. Chakravarty, J. L. Richardson, and 

D. J. Scalapino, Phys. Rev. B 24, 3682 (1981). 
[13] W. Janke, Phys. Rev. B 55, 3580 (1997). 
[14] E. Luijten and H. W. J. Blote, Int. J. Mod. Phys. C 6, 

359 (1995). 

[15] E. Luijten and H. Mefiingfeld, to be published; 
H. Mefiingfeld, M.Sc. Thesis, Johannes Gutenberg- 
University Mainz (unpublished). 

[16] D. R. Nelson and J. M. Kosterlitz, Phys. Rev. Lett. 39, 
1201 (1977). 

[17] D. J. Amit, Y. Y. Goldschmidt, and G. Grinstein, J. 

Phys. A 13, 585 (1980). 
[18] G. V. Matvienko, Theor. Math. Phys. 63, 635 (1985). 
[19] Z. Glumac and K. Uzelac, J. Phys. A 26, 5267 (1993). 
[20] M. Krech and E. Luijten, Phys. Rev. E 61, 2058 (2000). 




0.5 0.6 
K 

FIG. 1. The scaled order parameter ^ = K(m 2 ) of the in- 
verse-square Ising chain, as a function of the coupling K. The 
inset shows the estimates for lim_L_ 00 ty(K,L) for K > K c , 
along with a fit of the predicted singularity. 
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FIG. 2. The reduced magnetic susceptibility x = L d (m 2 ) 
of the inverse-square Ising chain as a function of the cou- 
pling, along with a (truncated) high-temperature series ex- 
pansion and a fit to the predicted asymptotic behavior 
e X p\Bo(K - K c )~ p + B 2 (K - K c f}. 
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FIG. 3. The critical correlation function g(r) of the in- 
verse-square Ising chain, for 5 < r < 5 ■ 10 s . The logarithmic 
decay can be followed for at least three decades in r. 
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FIG. 4. The decay exponent 9 of the spin-spin correlation 
function g(r) of the inverse-square Ising chain, as a function 
of the coupling. 
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FIG. 5. The specific heat of the inverse- square Ising chain, 
for a range of system sizes. Note that the maximum reaches 
its thermodynamic limit nonmonotonically. 



